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ABSTRACT 

We construct simple and useful approximation for the relativistic gas of massive 
i/-)" 1 particles. The equation of state is given by an elementary function and admits 

analytic solution of the Friedmann equation, including more complex cases when 
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the relativistic gas of massive particles is considered together with radiation or 



cr 



X 



with dominating cosmological constant. The model of relativistic gas may be 
interesting for the description of primordial Universe, especially as a candidate 
for the role of a Dark Matter. 
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During many years most of the works about cosmological solutions assumed stationary 
equation of state for the matter or vacuum sources of the Friedmann equations. The con- 
ventional equation of state is a linear relation between pressure and energy density P = wp. 
Different values of the parameter w correspond to different kinds of sources. For example, 
w = — 1 holds for the cosmological constant (vacuum energy), w = for the pressure- 
less (dust-like) matter, w = 1/3 for the radiation. Recently, due to the growing amount 
and quality of the observational/experimental cosmological data, theoretical considerations 
involved more complicated equations of state, with w depending on time or/and energy 
density. Perhaps the first examples of this kind were related to inflation, where the variable 
vacuum energy density was introduced due to the electroweak phase transition [I] and the 
non-trivial inflaton potential 

Theoretical investigations of the models with variable vacuum energy were fueled recently 
by the new precise measurements of the expansion rate of the Universe from the type la su- 
pernovae experiments jHj and also from the cosmic microwave background radiation £|] . The 
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existing data indicate that the Universe is mainly composed by the non-luminous sources, 
such as Dark Matter, responsible for 20-30% of the overall energy density balance, and the 
Dark Energy, responsible for 65-75%. One of the main candidates to the role of the most of 
the Dark Matter is a gas of weakly interacting massive particles, e.g., the ones corresponding 
to the broken supersymmetry. Let us notice that recent theoretical and phenomenological 
considerations of the supersymmetric neutralino do not rule out the light DM options |H] 
(see 0|B] for the review). The main candidate to be Dark Energy is the cosmological con- 
stant with equation of state w — — 1. The anthropic considerations show that the total value 
of the cosmological constant should be positive and in fact close to the observed one [3 EI] • 

However, since nobody can guarantee that the vacuum energy is red-shift independent, 
it is quite natural to meet a variety of alternative models for the Dark Energy, such as 
quintessence [H], Chaplygin gas |12j and the low-energy renormalization group running 
of the proper cosmological constant Most of these models lead to a variable w in 

the equation of state for the vacuum energy. This effect is achieved either by postulating 
this equation of state in the case of Chaplygin gas or by postulating a properly chosen 
quintessence potential. 

There is another possibility to meet the equation of state with a variable w, depending 
on the energy density and therefore on the red-shift. In the present article we shall consider 
the Universe filled by the ideal gas of relativistic massive particles 4 . This model may have 
interesting applications, e.g., in the early radiation-dominated Universe we can consider 
relativistic gas of massive particles as a model for the hot matter content. Furthermore, 
relativistic effects may be, in principle, relevant for the Dark Matter problem. Since we do 
not know exactly from what the Dark Matter is done, any possibility here deserves careful 
exploration. Indeed, when assuming that the Dark Matter is a relativistic gas, we suppose 
that it is composed (or has been composed in the earlier epochs) from the relatively light, 
massive particles weakly interacting with the baryonic matter and radiation. 

The equation of state for the relativistic ideal gas is known for a long time (see 
also, e.g., |IE]) 5 - The relation between pressure and energy density involves modified Bessel 
functions. Obviously, this form of equation of state in not very useful for cosmological 
applications. At the same time one can considerably simplify the cosmological model with 
relativistic gas without losing much of the physical sense. In order to do so we shall assume 
that, instead of following Maxwell distribution, all particles have equal kinetic energy. Below 
we call the model of relativistic gas of massive particles with equal energies the reduced 
relativistic gas. 

The "defect" in the equation of state which follows from the assumption of equal energies 

4 The cosmological model taking into account the relativistic effects related to the peculiar velocities of 
galaxies has been developed in [T3] , 

5 The relativistic version of Maxwell distribution follows from the corresponding generalization of the 
Boltzmann ii-theorem (see, e.g. |17|). 
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is not very significant. The numerical comparison with the Maxwell distribution is presented 
in the Appendix. It turns out that the difference between the two distributions does not 
exceed 2.5% even in the low-energy region, being negligible for the ultrarelativistic gas. Let 
us remember that the Maxwell distribution is also just an approximation to the real situation. 
For example, when considering the Maxwell distribution for the identical massive particles 
in the Early Universe we are disregarding interactions between these particles and radiation, 
and also differences between the masses of different kinds of particles. 

At the same time, the model of reduced relativistic gas provides a great advantage for 
cosmological applications. Starting from the reduced equation of state one can integrate the 
Friedmann equation analytically, leading to a nice and simple cosmological model interpo- 
lating between radiation-dominated and matter-dominated epochs of the Universe. In the 
present paper we shall develop this model and also consider more complicated cases of the 
Universe filled by reduced relativistic gas plus radiation and of the Universe dominated by 
the cosmological constant where matter content is the reduced relativistic gas. 

The paper is organized as follows. In the next Section we shall briefly describe the 
conventional model of relativistic gas and our model of the reduced relativistic gas. In Section 
3 we use the equation of state for the reduced model and obtain the scale dependence for 
the energy density and pressure. Section 4 is devoted to the solution of the cosmological 
model with the reduced relativistic gas for several interesting particular cases. In section 4 
we draw our conclusions. 



2 Reduced model for relativistic gas 

Consider a single relativistic particle with the rest mass m in a volume V. The dispersion 
relation for this particle has standard form 

e 2 - c 2 p 2 = m 2 c 4 , where p = mV = . (1) 

a/1 — v 2 /c 2 

An elementary consideration shows that the time average of the pressure produced by the 
particle on the walls of the vessel is 

~ 1 ^ (2) 



3U ./I - v 2 /c 2 ' 

For the gas of N such particles with equal kinetic energies e, we arrive at the following 
equation of state 

-2 N 2" 



P= P - 



where p = (3) 



is the energy density. Let us notice that w = P/p tends to 1/3 in the ultra- relativistic 
limit e — > oo and to zero in the non-relativistic limit e — > mc 2 . It proves useful to introduce 
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the following new notations: the density of the rest energy of the particles at the initial 
moment p\ = Nine 2 /Vq, where Vq is some fixed initial value of the volume, and also 
Pd = Pdiy) = PiVo/V = Nmc 2 /V, which shows how the same density evolves with the 
change of the volume. Using these notations we can cast the equation Q into the form 



P= P - 



l ~7 



(4) 



In order to understand better the difference between these formulas and the "correct" 
ones, let us consider the Maxwell distribution for the ideal gas of massive particles. The 
statistical integral for a single particle is given by the expression 

Z = J e- e / kT d 3 pd 3 q = 47tm 2 cV ■ ^(^) , (5) 

where K v (x) is a modified Bessel function of index v. The equation of state for the gas of 
N particles can be derived in a standard way 

pv = kTN il£v) T = NkT ' < 6 > 

while the average energy of the particle is 

The energy density p = Ne/V and pressure P are related by an implicit functional depen- 
dence © and (J7|), which has to be compared with the formula (jU for the reduced gas case. 
This comparison will be performed numerically in the Appendix, where we show that (@J) is 
an excellent approximation to the relations (JHJ), (J7J). For a while we conclude this section by 
an obvious observation that (jlj) is much simpler than (JHJ) plus (J2J). 

3 Scale dependence in the reduced model 

Let us use the equation of state (jlj) for the reduced relativistic gas and find how the 
energy density depends on the volume under adiabatic expansion. For this end we replace 
the equation (0J) into the conservation law 

dV dp 

In the cosmological setting the last equation implies that the reduced relativistic gas does not 
exchange energy with other entities like with radiation in the Early Universe or with baryonic 
matter and vacuum energy in the case when the reduced relativistic gas is considered as a 
model for the Dark Matter. 
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The differential equation 
P{V) 

where the initial condition has been defined as 

p(Vo) = [p\ + Pi] 



can be easily solved. The solution has the form 



V 

Pi \ y 



2 (V y/z 



1/2 



1/2 



(9) 



Remember that p\ = Nmc 2 /Vo is the rest energy density at the initial point V = Vq. 
Hence the second component p 2 can be interpreted as the energy density of the radiation 
component of the reduced relativistic gas. However the total energy density is not a simple 
sum of the two components but the square root of the sum of their squares. The non- 
relativistic or ultra-relativistic limits are achieved when one takes, correspondingly, p 2 = 
or pi = 0. It is easy to see that the expression © provides correct scaling laws in these two 
cases. 

For the sake of cosmological applications it is better to express the energy density as a 
function of the conformal factor a, where (a/ao) 3 = V/Vq. Then we arrive at the formula 

8 l 1/2 

) 



p(a) 



P\ 



OoV 

a ) 



+ pl( 



aoy 
a 



(10) 



with the initial condition p(a ) = [pf + p^\^ 2 - The density of the rest mass behaves like 
pd{a) = Pi(ao/a) 3 - The scale dependence of the pressure and parameter w are given by the 
eq. (HJ) and the equation 



w 



1 

3 



Pd( a ) 



(11) 

P'W J 

It is easy to see that the relations (jH) and (fTTj) predict dust-like (pressureless) scaling in the 
limit a — > oo and the radiation-like scaling w ~ 1/3 in the limit a —>■ 0. Hence our model 
of reduced relativistic gas can be regarded as an interpolation model between the radiation- 
dominated and matter-dominated evolutions. If the Universe, from the very beginning, were 
filled by ideal hot gases of massive particles, there would not be nuclear reactions and the 
scale behaviour of the energy density would be close to (fTUj). 



4 Solving the Friedmann equation 

Consider the cosmological model of the Universe filled by the reduced relativistic gas 
()10|). For the sake of generality, let us start by formulating the Friedmann equation for an 
arbitrary k and also include vacuum energy density, p\ = A/8tcG , and the radiation energy 
density p r (a) = p r0 /a 4 - In what follows we set ao = 1. The equation of interest has the 
form 

/d \ 2 k 8nG r , , . . -, „ . 

{-) +^ = -g-Ma)+PA + Pr(a)], (12) 
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where p(a) is given by (fTUj). One can easily add other matter sources and solve the resulting 
equation following the examples which will be elaborated below. 



4.1 Pure reduced gas model 

As a first particular case we shall integrate eq. (|T2j) for the pure reduced relativistic gas 
model, k = p\ = p r = 0. After introducing a new variable x = a 2 the equation becomes 

x 2 = — \fx + b , where b = ^ . (13) 

3 p{ 

The solution of this equation has the form 

.2 \ 3/4 



a 2 + ^|) =^/6nG Pl -t. (14) 
Pi J 

The last expression shows, again, that our model interpolates between the usual matter- 
dominated and radiation-dominated FRW solutions. In the non-relativistic case, p 2 = 0, we 
obtain directly the standard behaviour a(t) ~ t 2 ^ 3 . The ultrarelativistic regime can not be 
considered in a direct way, because the limit p 1 = is singular. Let us assume p\ <C p 2 
and expand until the lowest nontrivial order in the ratio (pi/ p 2 ) 



P2\3/2 

pJ 



4 \p 2 ) 



^/6nG Pl t . (15) 



The last expression is nothing else but the usual radiation-dominated solution with the 
shifted time variable 

a(t) = f^^El) ' Vt^To, where t Q = (16) 
\ o / V67rG • pi 

Indeed, such time shifts are irrelevant and we will not pay attention to them in what follows. 

In order to illustrate the behaviour of the conformal factor ()14[) in the transition period, 
let us choose the logarithmic parametrization for both time and conformal factor. For the 
sake of simplicity we fix t according to (|16j) and rewrite (|14j) in the form 



In 



e 2s , 

2 J 



pV \pJ 



(17) 



where e v = \JQixGp\ t and a = e s . 

The plot of v versus s depends on the value of a single parameter p 2 /pi and in all 
cases it clearly demonstrates the transition between the two linear asymptotic regimes. We 
present an example for p 2 j p\ = 1 at the Figure 1 . The reader can easily achieve similar 
plots for other values of this parameter, also in the presence of a radiation content. The 
solution for the last case is presented in the next subsection. 
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Figure 1: The log-log plot for the transition between the radiation-dominated 
and matter-dominated epochs in the pure reduced relativistic gas model. 



4.2 Radiation dominated epoch 



Consider the reduced relativistic massive gas in the radiation-dominated epoch. According 
to the known estimates, cosmological constant and space curvature are not very relevant in 
this case [TBI EH] and we can safely set k = p\ = 0. 

The solution of the Friedmann equation (|T2")) can be written as (parameter b is defined 
in (HI) 



4 
3 



Pi J 



3/2 



4p, 



rO 



Pi 



Pi . 



1/2 



32nGp i 



1/2 



t. 



(18) 



The last relation represents exact solution, but it is too complicated for the qualitative 
analysis. Let us consider the special situation when the radiation energy density is strongly 
dominated and the effect of the reduced relativistic massive gas is a small correction to the 
i 1 / 2 law. Our purpose is to evaluate this correction. For this end one has to expand the 
expression (here 7 = 1/2 or 3/2) 



pi J 



PrO 
Pi 



1 + ^V^Tb 

PrO 



(19) 



until the third order in the small parameter \/ a 2 + b . 



The result has the form 



i te) ^ + ^ ■ ( 



32iiGp r o 



1/2 



t. 



(20) 



where we disregard the initial value t . It is easy to see that the effect of the reduced 
relativistic gas is to accelerate the expansion of the Universe compared to the pure radiation 
content. The effect of relativistic gas of massive particles is weaker that the one coming 
from the dust-like matter with the same energy density. The illustrative plots of the reduced 
relativistic gas compared to radiation and dust cases are presented at the Figure 2. 
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Figure 2: Tie plot for the conformal factor a(t) in the following three cases: 
pure radiation a(t) ~ t x l 2 is represented by the dashed line, pure dust-like case 
by the dots line, the case of pure reduced relativistic gas is represented by the 
continuous line. 



4.3 Cosmological constant dominated epoch 

The next relevant particular case is the reduced relativistic gas as a model for the Dark 
Matter, in a Universe dominated by the vacuum energy. After the usual change of variable, 
x = a 2 , we arrive at the solution in the form of an integral 



dx 32tvG , . 

-5-t, (21) 



[p\x 2 + pi\?x + b\ 

with b defined in (|13|). This integral is rather complicated and difficult to evaluate ana- 
lytically. Hence we can either apply the numerical method or use the dominant role of the 
vacuum energy density. The results of numerical analysis for = 0,7 and Qm = 0,3 are 
shown at the Figure 3. The presence of the reduced relativistic massive gas results in the 
slower expansion of the Universe compared to the pure cosmological constant case. 

Consider an approximate analytical solution based on the assumption of the dominant 
contribution of the vacuum energy density. We rewrite the integral in the l.h.s. of eq. ()21j) 
in the form 



1 f dx pi \Jx + b 

where 7 = — 



y/pX J x a/1 + 7 Pa x<1 

and expand into power series in the small parameter 7. Then the integration becomes trivial 
and taking the lowest nontrivial order into account we arrive at the solution 

In a + X[a) = At, (22) 

where 



y/^F+b y/^Tb 1 , /y/d^Tb-Vb 
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Figure 3: The plot of conformal factor a(t) for the following three cases: cos- 
mological constant plus ordinary pressureless matter (points line), cosmological 
constant plus pure radiation (dashed line) and cosmological constant plus reduced 
relativistic gas (continuous line). One can see that the effect of the reduced rel- 
ativistic gas is intermediate between the ones of the dust-like and radiation-like 
matter contents. At later times all three curves converge and only the cosmolog- 
ical constant remains relevant. 



is a small term and A = a/87tGpa/3. One can easily find an approximate explicit formula 
for a(t) starting from the expression 

a{t) = [1 + /(*)] • e x \ |/(f)l<l, (24) 

and using the smallness of X[a\. Replacing (|24|) into (|22|) we arrive at the solution 
/(£) = — X[e A *]. This expression can be replaced into (|2*4*jl to give 

a(t) = e xt (l-X[e A ']) . (25) 

Qualitatively the behaviour of a(t) fits with the numerical analysis. The presence of reduced 
relativistic gas slows the acceleration of the Universe caused by p A . 

5 Conclusions 

We constructed a simple and useful model of reduced relativistic gas of massive particles, 
starting from the "primitive distribution" for kinetic energies. Our model is a very good 
approximation to the much more complicated Maxwell distribution (see Appendix). The 
main advantage of the reduced model is that it admits analytical derivation of the dependence 
p(a) in terms of elementary (and very simple) functions. Moreover, one can easily integrate 
the Friedmann equation for the Universe filled by the reduced relativistic gas, again the 
solution is given by elementary function. In this way we arrive at the interpolation between 
the cosmological solutions for the matter-dominated and radiation-dominated cases. 
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In the more complex situation, when the sources of the Friedmann equation include also 
radiation or the cosmological constant term, one can obtain the solution in the form of an 
integral. The integration can be performed numerically or analytically using the assumption 
of radiation dominance or cosmological constant dominance. In the last case the reduced 
relativistic gas can be viewed as a model for the Dark Matter. 

The model of reduced relativistic gas may be, in principle, testable. Imagine the Dark 
Matter is composed by the relatively light, weakly interacted massive particles which have, 
at present, w ~ equation of state. One can not rule out the possibility that the relativistic 
effects of these particles were relevant in the earlier epochs of the Universe, e.g. in the struc- 
ture formation period. Then, some the traces of these relativistic effects may be eventually 
found in the precise CMB measurements. Therefore it would be interesting to investigate 
density and metric perturbations in this model. We postpone this problem for the future 
work. 



Appendix 

The purpose of this Appendix is to compare the results (jOJ and ((Zj) of the Maxwell distribution 
for the relativistic gas of massive particles and the corresponding relation for the reduced 
relativistic gas model. In the formulas © and ((7j) the temperature kT plays the role of 
parameter and we are in fact interested only in the dependence between P and p. Hence 
we solve (JEJ) with respect to pressure and replace it into ((Zj) 

, P N K 3(pd/P) p , 0R x 

^ p) = KMPj^- p ' (26) 

where the subscript M indicated Maxwell distribution. 

It is easy to see that the non-relativistic limit P — > and the ultrarelativistic limit 
P — > p/3 of the last relation coincide with the one for the reduced model (jlj). In order to 
compare the two expressions numerically at the intermediate scales, let us rewrite (@J) in the 
form similar to (J26|) 



P(P) = \P+ y^ + Pl (27) 

After assigning numerical value to pi (we set p\ = 1, but the result does not depend on this 
choice.), one can plot the relative deviation 



8p 



\P~ Pm\ 



Pm 

versus P. The plot for this dependence in presented at the Figure 4. It is easy to see that 
the relative deviation S p achieves its maximum of about 2.5% in the non-relativistic region 
and becomes completely negligible in the higher energy region. 
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